RATIONAL RUIJSENAARS-SCHNEIDER HIERARCHY 
AND BISPECTRAL DIFFERENCE OPERATORS 



PLAMEN ILIEV 

Abstract. We show that a monic polynomial in a discrete variable n, 
with coefficients depending on time variables ti,t2, ■ ■ ■ is a r-function for 
the discrete Kadomtsev-Petviashvili hierarchy if and only if the motion 
of its zeros is governed by a hierarchy of Ruijsenaars-Schneider systems. 
These r-functions were considered in [TS], where it was proved that they 
parametrize rank one solutions to a difference-differential version of the 
bispectral problem. 



1. Introduction 

In [4J, Airault, McKean and Moser discovered a mysterious connection be- 
tween equations of KdV type and the Calogero- Moser system. They showed 
that the motion of the poles of a rational solution to the KdV or Boussi- 
nesq equation that vanishes at infinity is described by the Calogero-Moser 
system |6J, with some constraint on the configuration of poles. Krichever 
[16| observed that the poles of the rational solutions to the KP equation 
that vanish at x = oo, move according to the Calogero-Moser system with 
no constraint. Shiota [22] extended this phenomenon to the whole KP hi- 
erarchy, which combined with the work of Adler [Ij led to a simple explicit 
formula for the r-function. 

A surprising link to the above theory was observed by Duistermaat and 
Griinbaum ^ in connection with a problem in limited angle tomography 
|10j . known as the bispectral problem. As originally formulated, this prob- 
lem asks for which ordinary differential operators L{x,d/dx) there exists 
a family of eigenfunctions ^(x,z) that are also eigenfunctions for another 
differential operator B{z,d/dz) in the "spectral parameter" z. In the case 
when the operator L{x,d/dx) belongs to a rank one commutative ring of 
differential operators (i.e. L commutes with an operator of odd order), the 
solution of the bispectral problem (up to translations and rescalings of x and 
z) are precisely the operators which can be obtained by finitely many ratio- 
nal Darboux transformations from Lq = d^/dx^. This combined with work 
of Adler and Moser p] shows that the rank one solutions of the bispectral 
problem are exactly the rational solutions discovered in [4j . Wilson [26j pro- 
posed to extend the problem to commutative rings of differential operators. 
Such a ring is called bispectral when there is a joint eigenfunction of the op- 
erators in the ring that is also a joint eigenfunction of a ring of differential 
operators in the spectral variable. An important invariant of such a ring is 
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its rank, meaning the dimension of the common space of eigenfunctions to 
the operators belonging to the ring. He proved that the bispectral maximal 
rank one commutative rings of differential operators are parametrized by a 
sub-Grassmannian Gr^'^ of Sato's Grassmannian [2T], which corresponds to 
the rational solutions of the KP equation studied by Krichever [16] . More- 
over, in a subsequent paper [27], Wilson gave a beautiful explanation of 
the bispectral property based on the connection with Calogero-Moser sys- 
tems and their geometric description [l5j. He also deepened the mystery 
by showing that the correspondence between the Calogero-Moser and the 
KP systems extends even to the locus where the particles collide with each 
other. For a very nice characterization of the Grassmannian Gr^'^ in terms of 
representation theory see the recent work of Horozov [13]. For an intriguing 
connection of the above theory to noncommutative geometry see [5]. 

In [12], jointly with Luc Haine, we constructed rank one commutative 
rings of difference operators in a discrete variable n G Z, corresponding to 
a flag of nested subspaces, each of which belongs to Gr'^'^. We showed that 
the common eigenfunction of the operators in the ring is also the common 
eigenfunction of a maximal rank one commutative ring of differential op- 
erators in the spectral variable, i.e. they provide rank one solutions to a 
difference-differential version of the bispectral problem. The corresponding 
r-functions r(n; t) are polynomials in n and give rational solutions of the 
discrete KP hierarchy. 

In the present paper, we investigate the motion of the zeros of polynomial 
(in n) r-functions of the discrete KP hierarchy. We show that a monic 
polynomial in n of degree is a r-function for the discrete KP hierarchy 
if and only if the motion of its roots {xi}^^ is governed by a hierarchy 
of Ruijsenaars-Schneider systems. We restrict out attention to the generic 
situation when the roots satisfy the constraints Xi — Xj ^ {0, 1} for i ^ 
j. This condition means that the rational solution to the simplest zero 
curvature (Zakharov-Shabat) equation for the discrete KP hierarchy has 2N 
distinct poles. It is a challenging problem to investigate the more general 
case allowing collisions of the poles. 

The paper can be thought of as a discrete analog of Shiota's paper [22] . 
In particular, from the proof, we can easily write an explicit formula for the 
r-functions in terms of the Ruijsenaars-Schneider hierarchy, which implies 
that they parametrize the rank one solutions of the difference-differential 
version of the bispectral problem constructed in [12] . 

We note that there is a related work of van Diejen [25] in the case of 
second-order difference operators, where the dynamics of the zeros of the 
solitonic Baker-Akhiezer function in the spectral variable z is studied. For 
soliton solutions of KP and 2D Toda equations, see Ruijsenaars [19] and 
van Diejen-Puschmann [25] . and for elliptic generalizations see Krichever- 
Zabrodin [T7]. For a g-deformation of the KP hierarchy and connections 
with the bispectral problem see [14] . 

The paper is organized as follows. In the next section, we briefly introduce 
the necessary ingredients of the discrete KP hierarchy. The approach follows 
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closely [12] . which leads to r-functions that differ by an exponential factor 
from the ones constructed in [H [23] . In Section [3] we formulate the main 
result of the paper and its connection to the bispectral problem. For a very 
nice account on the difference-differential version of the bispectral problem 
and its relations to orthogonal polynomials and the Toda lattice see [TT]. 
Section m is devoted to the proof of the main result. 



2. The discrete KP hierarchy and t-function 

We denote by A and V the customary forward and backward difference 
operators acting on functions of a discrete variable n € Z by 

A/(n) = f{n + 1) - f{n) and V/(n) = /(n) - f{n - 1). 

The formal adjoint to A is A* = —V. If we define 

A^' • fin) = Yr.) i^'f){n + j - i)A^-\ for ah j G Z, 
i=o 

we obtain an associative ring of formal pseudo-difference operators 



R{A} = ix= ^ a,{n)AA. 

[_ j=—oo ) 



We denote by X+ = X^^=o o,j{n)A^ the positive difference part of X and by 
X_ = Ylij=-co 0'j{n)A^ , the Volterra part of X. 

The discrete Kadomtsev-Petviashvili hierarchy (in short KP) is the family 
of evolution equations in infinitely many time variables t = {ti,t2,ts, . . .) 
given by the Lax equations 

Sr^-Hn..L], (2.1) 

where L is a general formal pseudo-difference operator of the form 

oo 

L = A + ^aj(n)A"-''. 

j=0 

A r-function for the hierarchy (j2.ip can be defined as follows. First, we 
define a wave operator 

oo 

W{n;t) = 1 + Ywj{n;t)A-^ , 

which conjugates L to A, that is 

L = WAW-'^. (2.2) 
The vector fields (|2.1|) can be extended by 

dW I. 
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For simplicity, we denote by Exp(n; t, z) the exponential function 

Exp(n; t, z) = (1 + z f exp Uz"^ . 

The wave function w{n; t, z) and the adjoint wave function w*{n; t, z) of the 
discrete KP hierarchy (12. 1|) are defined by 

w{n] t, z) = W{n; t)Exp(n; t, z) 

1 + ^ ' ' + \' ' +■■■] Exp n; t, z) 2.4a 
z z'^ J 

and 

u;*(n;i,z) = {W~'^{n - 1; t))* Exp"^(n; t, z) 

1 + + + • • • ) Exp-^n; t, (2.4b) 

z z'^ J 

The functions w{n;t,z) and i(;*(n;t, z) can be written in terms of a r- 
function as follows 

w{n; t, z) = '^(^'^ -[^ ^'' Exp(n; t, z), (2.5a) 

and 

it;*(n; t, z) = I^!^ii±J^JlExp-i(n; t, z), (2.5b) 

where [z] = (z, 2:^/2, z^/3, . . . ). We refer the reader to [12] for more details 
and proofs of the above construction. 

Remark 2.1. It is well known (see for instance [71 Proposition 5.1.4, p. 75]) 
that the Lax equations (j2.ip imply the zero curvature (Zakharov-Shabat) 
equations 

^-^^ = [(^™)+,(^'=)+], (2.6) 

where k,m G N. In the differential case, the first flow ti corresponds to 
a translation in the spatial variable and therefore the simplest (nontrivial) 
zero curvature equation can be obtained for m = 2 and A: = 3. This leads 
to the KP equation, which gave the name of the whole hierarchy. In the 
discrete case, the first flow is no longer trivial and the simplest zero curvature 
equation will correspond to the choice m = 2 and k = 1. In the example 
below we carry out this computation explicitly, which leads to a nonlinear 
partial differential-difference equation for the function ao(n;t). 

Example 2.2. Let us take m = 2 and A; = 1 in (j2.6p . Clearly, (L)^ = 
A + ao(n;t) and a short computation shows that 

{L^)+ = + (ao(n; t) + ao(n + 1, t)) A + ag(n; t) + ai(n, t) + ai(n + 1, t). 
Next, we see that 

[{L^)+, (L)+] = {A^oin, t)-A (ai(n; t) + ai(n + 1; t))) (A + 1) . 
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Using the relations above and comparing the coefficients of A* for i = 0,1 
in (j2.6p with m = 2 and A: = 1 we get the system 

d (ao(n; t) + ao{n + 1; t)) 

= A^ao{n;t) - A {ai{n;t) +ai(n + l;t)) 
dao{n;t) d {a^in; t) + ai{n;t) + ai(n + l;t)) 

= A^ao(n;t) - A {ai{n;t) + ai{n + l;t)) . 
EUminating ai(n;t) we obtain the following equation for = ao(n;t) 

— Aao = ^ (Aag - 2Aao) + (Aoq + 2ao) . (2.7) 
012 oil oil 



3. Polynomial t-functions and the dynamics of their zeros 
The main result of the paper is the following theorem. 

Theorem 3.1. Let xi{t),X2{t), . . . ,XN{t) be smooth functions oft = (ti,t2i • • • 
such that Xiit) — Xj{t) ^ {0, 1} for i ^ j and dxi{t)/dti ^ in a neighbor- 
hood oft = 0. Let us define functions yi{t),y2{t), . . . ,yN{t) by the following 
relation 



e 



-y^{t) 



dxi{t) Xi{t) - Xs{t) 



dti fJ^Xi{t) - Xs{t) + 1 



Then the following conditions are equivalent. 
(i) The function 



N 



r{n;t)=ll{n-Xi{t)), (3.2) 



i=l 

is a T -function for the discrete KP hierarchy (j2.ip . 
(ii) The motion of {xi{t),yi{t)}j^^ is governed by the Ruijsenaars- Schneider 
hierarchy of Hamiltonian systems 

J^(;:)-<-')'(-T4'/t.). (-) 

where Hk = tr(y*^), and Y is an N x N matrix with entries 

-y ^ _ _L 1 

y^J = + ' ' , n ^ ■ (3-4) 

s=l 

Remark 3.2. Although the Hamiltonian system above differs slightly from 
the standard rational Ruijsenaars-Schneider system, one can easily connect 
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the two. Indeed, the rational Ruijsenaars-Schneider model is a dynamical 
system, whose equations of motion can be written in the following form 

N 2 

see |20l formulas (B22)-(B23), p. 402]. If we now put qj = i'yxj, then 
(jS.Sp gives precisely the dynamical system p.3p for the first flow d/dti, see 
equation (j4.7p . 

Remark 3.3. The solutions of the KP equation 

3 r 1. ,] 

-Uyy = I lit - -{UxxX + QUUx)j 

that are rational in x and vanish as x — > oo have the form 

N ^ 

When all Xj are distinct their motion is governed by the Calogero-Moser 
system, as shown in [16]. The case discussed in the present paper is similar in 
the following sense: the solutions described in Theorem 13.11 provide rational 
solutions of equation (12. 7p . which vanish as n ^ oo. These solutions have 
simple poles (as functions of n) at the points {xj,Xj — IjjLn see formula 
(j4.3p . The condition Xi — xj ^ {0, 1} for i ^ j simply means that all these 
poles are distinct. 

As a consequence of the proof of Theorem 13.11 we also obtain an explicit 
formula for the r-function in terms of {xi,yi}^i at ti = ^2 = • • • = 0. Let 
us denote by X the diagonal matrix with entries Xi{t), i.e. 

X = diag(xi(t),X2(t),... ,XN{t)). (3.6) 

Corollary 3.4. Let X^ and be the matrices X and Y, defined by (j3.6p 
and (|3.4p at ti = t2 = ■ ■ ■ = 0. Then the r-function in equation ()3.2p can 
be computed from the following formula 

oo 

r(n; t) = det (nl - X^ + Y,jW " Y^){-Yy~^) , (3.7) 

i=i 

where I is the identity N x N matrix. 

From formula (j3.7p it is easy to see that r(n; t) = r(0; ti + n,t2 — n/2,t3 + 
n/3, . . . ), where r(0; t) is a r-function for the (continuous) KP hierarchy, cor- 
responding to a plane in Wilson's adelic Grassmannian. Thus, the results 
in [12] imply that the functions r(n; t) described in Theorem 13.11 parame- 
trize rank-one solutions to a difference-differential version of the bispectral 
problem. More precisely, there exist a rank-one commutative ring A of dif- 
ference operators in the variable n, and a rank-one commutative ring A' of 



RUIJSENAARS-SCHNEIDER HIERARCHY AND BISPECTRAL OPERATORS 7 



differential operators in such that 

Cw{n; t, z) = fc{z)w{n; t, z), V£ G A 
Bw{n;t, z) = gts{n)'w{n;t, z), £ A' 

where fc{z) and gein) are functions of z and n, respectively, and w{n;t,z) 
is the wave function defined by (|2.5ap . 



4. Proof of Theorem 3.1 

The strategy of the proof is as follows. For the implication (i)=^(ii), we 
investigate the motion of the poles with respect to the first fiow d/dti and 
we write the corresponding dynamical system in an appropriate Lax form. 
This represents a discrete analog of some of the results in [T8] , except 
that in the continuous case the first nontrivial flow is d/dt2. Next, we adapt 
the approach in [22j to establish the Hamiltonian equations for the higher 
flows d/dtk, k > 2. The opposite direction can be deduced by using the 
connection between polynomial (in ti) r-functions of the KP hierarchy and 
polynomial (in n) r-functions of the discrete KP hierarchy [12]. 

Let us start with the implication (i)=^(ii). Prom equations (j2.4p . (|2.5p 
and (|3.2p it is clear that we can write Wk{n;t) and ■w'^{n;t) as 



N 

y 

^ n - Xi{t) 
1=1 



In particular, for /c = 1 we see that wi,j(t) = dxi{t)/dti and w\^{t) = 
-dxi{t)/dti. From ([22]), ([O]) and ([Tla]) it follows that 

^^^^^Af) = (A + ao(n; t))w{n- 1, z). (4.2) 

Writing (12. 2p as LW = W IS. and comparing the coefficients of A*^ on both 
sides gives 



ao(n; t) = —wi{n + l;t) + wi{n; t) 

N 



E 



dxi{t) 



{n - Xi{t)){n + I - Xi{t)) dti 



(4.3) 
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where in the last equahty we used ()4.1ap for k = 1. Plugging the last formula 
for ao(n;t) in (I4.2jl and using (j4.1ap we get the following identity 

/ Wk+i,iit) _^ 1 dwk,i{t) ^ Wk,i{t) dxi{t) \ 
^-^\n — Xi{t) n — Xi{t) dti {n — Xi{t)Y dti j 

I Wk+l,i{t) Wk,i{t) Wk,iit) \ 

^^yn + l-Xiit) n + l-x^{t) n-Xi{t)J 

( I dxi{t) \ ( Wk^i{t) \ 

y2^^{n-Xi{t)){n + l-x^{t)) dh ]\^^n-Xi{t)j' 

Notice that (14.41) can be rewritten as a polynomial identity in n, which is true 
for every n G Z and therefore, it will be true for every n G C. Computing 
the residue at n = Xi{t) — 1 we obtain 

^ Wkjjt) dxijt) 

If we denote Wkit) = {wk,i{t),Wk,2{t), ■ ■ ■ ,Wk,N{t)Y , e = (1,1,...,!)*, then 
the last formula can be rewritten in vector notations as Wk+i{t) = {—Y)wf;{t), 
where Y is the matrix defined in Theorem 13.11 Thus we see that 

w,{t) = (-y)'=-i|-e, (4.6) 

where X is the diagonal matrix given in equation (j3.6p . 

Computing also the residue of equation (j4.4p at n = Xi{t) we obtain 

/ dxi{t) ^ 1 dxj{t)\ 

Using the last identity and (j4.5p we can eliminate Wk+i,i{t). For k = 1 this 
leads to the following second-order differential equation for Xi{t) 

d'^xiit) _ _^dxi{t) 



dti dti 

N 

^ - Xj{t)){x,{t) - Xj{t) + l){Xi{t) - Xj{t) - 1) dti 



^ 1 dxj{t) (4-7) 



which will be needed later. 

Similarly, if we work with the adjoint wave function w*{n;t,z) we can 
show that it satisfies the following equation 

dw:(^^ = (V - ao(n - l;t))w*in; t, z). 
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Denoting wl{t) = {wl ^{t), wl 2(t), . . . ,wl. j\f{t)Y we obtain as above that 

wtit) = -^^{-rr-'e. (4.8) 
Next, we use (|2.3p . From equations (j2.2p and (|2.4p we deduce that 

oo 

L'' =W{n;t)A^W{n;t)-^ = ^ Wi{n;t)A''-'-^ ■ w*{n + l;t), 

i,j=0 

where wo{n;t) = tt;g(n;t) = 1. This shows that 
fc+i 

j=0 

On the other hand 

Win; t) = 1 + f f; ' ^Jl^'A A-^ + 0(A"^). 

Plugging the last two formulas in (j2.3p . and equating the coefficients of 
on both sides we get 

/ 1 d^Xijt) 1 dxijt) dx,{t) \ 

^^yn-Xiit) dhdtk {n-Xi{t)y dti dtk J 
fc+i 

= - Wk+i-j{n; t)w*{n; t). 

j=0 

The last equality holds for every n G Z and therefore it must hold for every 
n G C. Comparing the coefficients of (n — Xi(t))~'^ gives 

dxjjt) dxjjt) A ^ 

Let us denote by li the elementary N x N matrix having 1 at entry (i, i) 
and everywhere else. Using the last identity, (j4.6p and (j4.8p we obtain 

dxAt) dxAt) / X ^ ^ 

i=i 

k 

Notice that dX/dtili = dxi{t) / dtili and therefore, we can cancel dxi{t) / dti 
and the last formula reduces to 

OTfc ^ oti 
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For every NxN matrix A we have e^Ae = tr(Aee*), and thus we can rewrite 
the last formula for dxi (t) / dt^ as follows 



^ = (-i)^-Hr( j:y-^/.y-^&e 

dtk \ ^ dti 

\j=i 



From the definitions of matrices X and Y it is easy to see that ^ee* = 
— {XY — YX — Y + I). Making this substitution and using the fact that 
tr{AB) = ir{BA) we get 

dXi{t) __ /^yi-lj^yfc-i(^y _YX-Y + I)] 



dU 



0=1 



= (-l)'=tr ^ (^iY^-^XY^ - liY'^-^+^XY^-^ - hY^ + hY^^^ 

= {-lfiv{li{XY'' - Y^X - kY^ + kY^-^)). 

However, UX = Xli and therefore tr(/jXy*=) = tr(/jy*^X). This shows 
that 

^ = k{-lf tr((/, - hY)Y^-'). (4.9) 

On the other hand it is easy to see that 

dY 



„ h- liY. 

oyi 



Thus 



d_ 

dy-, 



= ktrl^Y'-A =ktT{{h-hY)Y^'^). 



The last formula combined with (j4.9p gives the first equation in (j3.3p . In 
order to prove that the second equation holds, we first notice that (j4.7p is 
equivalent to the Lax equation 

§^ = [^,M], (4.10) 
where M is an x matrix with entries 

M-ij = -— — for i 7^ 7 

Xi{t)-xj{t) dti 

^ _ 1 dxkjt) _ 1 dxk{t) 

f^^Xi{t)-Xk{t) + l dti f^^Xiit) -Xk{t) dti 
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Differentiating ()3.ip with respect to tk we get 



dyi{t) _ I (^^ii*) 1 , "9 



dtk dtk 



dti 



, s 



Xi{t) - Xsjt) + 1 
L Xi{t) - Xs{t) 



(4.11) 



The derivative in the second term on the right-hand side of (j4.1ip can be 
evaluated using (|4.9p 



d 



Xjjt) - Xsjt) + 1 

L Xi{t)-Xsit) 



N 



N 



\I,{I-Y)Y 



k-l 



^ s=l 



For the first term we use both (|4.9p and (j4.10p : 



d 



log 



dxi{t) 



dtk \ dti 



k{-l) 



k(dx,{t)^-i d 



dti J dt] 



dxjjt) 
dti 



k/dxi{t) 



^ I k~\ k N 



dti 



tr {Mh{I -¥)-{!- Y)I,M) Y 



fc-i 



To simplify the formulas, let us denote Y = I — Y and {M)ij = {l — Sij)Mij 
(i.e. M is the matrix obtained from M by replacing the diagonal entries with 
zeros). Then, the last two formulas combined with (|4.1ip show that 

%(t) 



k{-l)''+hiiBY 



dtk 



k-U 



(4.12) 



where 



B 



dxjjt) 
dti 



-1 



N 



N 



[MUY - YUM) - E ( ^ log n 

J = l \ -I S = l 



lIjY. 



A straightforward computation now shows that 

N , N 



B 



dY 

dxi 



E — —h-y:— 



-I,,Y 



which combined with ()4.12p gives 



%(0 
dtk 



ay. 



A;(-l)'=+^tr —Y 
\ dxi 



k-l 



-i)'+'^tr(y*=), 
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completing the proof of (13. 3p . 

Conversely, assume now that (IS.Sp holds. Let us consider {xi{t),yi{t)}'j^i 
and the corresponding matrices X and Y at the initial time ti = t2 = • ■ ■ = 
and let us denote x° = Xj(0), y° = yj(0), X° = X|t=o, = y\t=o- Notice 
that 

rank(x°y° - + / - y°) = 1. 

Using the Cauchy determinant formula we see that 

det(/-y°) = e~^^=^y^ / 0. 

Thus, if we denote = X^{I — Y^)^^ we have 

rank([X°,y°] +/) = 1. 

Therefore the pair (X^, Y^) defines a plane in Wilson's adelic Grassmannian 
Gr*^*^, see [27] , The corresponding r- function can be computed by Shiota's 
formula 

oo 

fO(t) = det ( - + J^jt,(-yOy-i), 
i=i 

see [221 Corollary 1, p. 5845]. Applying [HI Theorem 2.4, p. 290] we deduce 
that f^(ti + n,t2 — n/2,t3 + n/3, . . . ) is a r-function for the discrete KP 
hierarchy (j2.ip . Multiplying by the nonzero constant factor det(/ — Y^) we 
see that 

f (n; t) = det(/ - y°)f ° (ti + n, is - ^, ^3 + ^, • • • ) 

oo 

= det (ni -x^ + ^ - y°)(-y°)^-^) 

is a r-function for the discrete KP hierarchy. Clearly, f{n;t) is a monic 
polynomial in n, and therefore, by the first part of theorem, its roots Xi{t) 
and the corresponding yi{t) will satisfy the Hamiltonian systems (13. Sp . To 
complete the proof we show that r(n; t) given by (j3.2p coincides with f(n; t) 
defined above. Since the roots of r(n; t) and f(n; satisfy the same systems 
(j3.3p . it is enough to show that x'j = Xj{0) and = yj{0). This follows 
easily from the explicit formula for f{n; t): 

f(n; ti, 0, . . . ) = det(n/ - X° + - y°)) 

= n - - + n ^^ozV + 

j=i \ s=i i / 

Remark 4.1. It would be interesting to see if one can use the explicit 
formulas for r-functions of KP hierarchy in terms of matrices satisfying rank 
one conditions [9j and the construction of r-functions for g-KP hierarchy 
from classical ones to extend the above proof and to show that every solution 
of the g-deformed Calogero-Moser hierarchy described in |14t Theorem 6.1] 
leads to a r-function for g-KP. This would give a one to one correspondence 
between rational solutions to the g-KP hierarchy (which also parametrize 
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rank one solutions to a bispectral problem for g-difference operators) and 
q'-deformed Calogero-Moser type systems. 

Acknowledgments. I thank a referee for his suggestions to improve an 
earlier version of the paper. 
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